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ABSTRACT 



We construct a mass model for the spiral lens galaxy 2237+0305, at redshift 
z;=0.04, based on gravitational-lensing constraints, HI rotation, and new stellar- 
kinematic information, based on data taken with the ESI spectrograph on the 10m 
Keck-II Telescope. High resolution rotation curves and velocity dispersion profiles 
along two perpendicular directions, close to the major and minor axes of the lens 
^3 galaxy, were obtained by fitting the Mgb-Fe absorption line region. The stellar rota- 

tion curve rises slowly and flattens at r ~ 1.5" (~1.1 kpc). The velocity dispersion 
£NJ ■ profile is approximately flat. A combination of photometric, kinematic and lensing 

' information is used to construct a mass model for the four major mass components of 

the system — the dark matter halo, disc, bulge, and bar. The best-fitting solution has 
a dark matter halo with a logarithmic inner density slope of 7=0.9±0.3 for pdm c< r 1 , 
a bulge with M/Lb=6.6±0.3T q , and a disc with M/Lb =1.2±0.3T©, in agreement 
with measurements of late-type spirals. The bulge dominates support in the inner 
£N| ■ regions where the multiple images are located and is therefore tightly constrained by 

the observations. The disc is sub-maximal and contributes 45 ± 11 per cent of the 
rotational support of the galaxy at 2.2r<i. The halo mass is (2.0 ±0.6) x 10 12 M Q , and 
the stellar to virial mass ratio is 7.0 ± 2.3 per cent, consistent with typical galaxies of 
the same mass. 
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of lensing information with photometric and kinematic 
data can be used to break degeneracies inherent to each 
source of information (e.g. the mass-sheet degeneracy 
and the mass-profile anisotropy degeneracy) to pro- 



Kooi 


Dmans et al. 1998; Winn et al. 


l2003l;lKooDmans & Treul 


2003 


; Treu & KooDmansI 2004 


; ICzoske et all 2008; 



in spiral galaxies which have multiple mass components 
that have to be modelled simultaneously, and the limi- 
tations of individual methods are particularly stringent. 
Specifically, in the absence of kinematic information, lens 
image positions can accurately constrain the enclosed mass, 
and to some extent also constrai n the mass distr i bution 
up to a constant mass-sheet dGorenstein et all 1 19881 : 
IWambsganss fe Paczvnskl 1 19941 ; IWucknitj |20Q2j ). When 
assuming a stellar mass-to-light ratio for the different opti- 
cal components (e.g. disc and bulge), based, for example, 
on stellar-population synthesis models, the dark matter 
contribution to the galaxy can then be inferred. In the 
absence of lensing information, but with kinematic data in 
hand, the overall galactic mass distribution is known, but 
again a stellar mass-to-light ratio must be assumed for the 
optical components in order to infer the contribution of the 
dark matter (the dis c-halo degeneracy, IPerez et al.l |2004| ; 
IVerheiien et al. 1 120071 ). Because of the rising contribution 
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of the stellar mass components, this becomes increasingly 
more difficult towards the inner regions of the galaxy. With 
the combination of lensing and high-resolution kinematic 
information, however, the stellar mass-to-light ratios of the 
optical components can be further constrained, breaking 
some of the de generacies between mass, mass-slope and 
aniso tropv (e.g. ITrott fc Websterl |2002| ; | Treu fc Koopmansl 
|2004| . and references therein). 

An accurate model for the mass distribution in spiral 
galaxies can provide interesting information about both the 
baryonic and the non-baryonic components. Disentangling 
these components provides information about the mass-to- 
light ratios as well as the influence, shape and density pro- 
file of the dark matter. The main motivations of this work 
are to: (1) break the classical degeneracy between the stel- 
lar component and the dark -matter halo in spiral galaxies 
l|van Albada fc Sancisi 1 9861 ). using the additional informa- 
tion that lensing provides; (2) determine independent mass- 
to-light ratios for the bulge and disc components; ( 3) deter- 
mine whether the disc compo nent is "maximal" (e.g. lSackettl 
Il997l : ICourteau fc~Rrxl Il999h . In more general terms, un- 
derstanding the relative mass of the different components 
can help constrain mode ls for disc galaxy formation (e.g. 



Dutton et al. II2OO5I, 120071 . and references therein, see also 



Shin fc Evansl ( 20071 )). In addition, the presence of massive 



discs in lens galaxies gives rise to interesting perturbations to 
the standard elliptical models adopted for early-type galax- 
ies with a number of interesti ng consequences for the topol- 
ogy of the time delay surface (|Moller et al.|[20o3 ). 

In this paper, we address these issues by studying the 
gravitational lens galaxy 2237+0305. Newly measured rota- 
tion curves and line-of-sight velocity dispersion profiles along 
two directions close to the major and minor axes - obtained 
from Keck-ESI spectroscopy - are combined with photomet- 
ric and lensing data from HST to construct a detailed mass 
model for the system. This is used to separate luminous and 
dark matter over scales ranging from the central parts dom- 
inated by the bulge out to several exponential scale lengths 
of the disc component, with the aid of HI data. 

The lens system 2237+0305 - chosen for this study due 
to the low redshift of the lens galaxy zi = 0.04 - has been ob- 
serve d extensively since its discovery in 1985 (|Huchra et al.l 
1 19851 ). In addition to providing i nformation about the 
global mass s t ructu r e of the gala x y dSchneider et al] 1988; 
iKent fc Falccl 1 19881 ; iMihovl l200ll ; ITrott fc Webster! |2002| ). 
it has also been used to study the properties of stars, or 
compact objects, in the lensing gal axy through microlens- 
ing (e.g. iGil-Merino fc Lewis! 12003 ) . and the structure of 
the broad- line region of t he lensed quasar dWayth et all 
120051 ; IVakulik et all 120071 ; lEigenbrod et~al] l2008bl ). PrevT 
ously published spectra of 2237+0305 have either concen- 
trated on studying the quasar spectrum, or only presented 
a small wavel ength range for the galaxy spectrum. The dis- 
covery paper. iHuchra eta l. ( 1985), presented a quasar spec- 
trum taken with the MMT and featured both the CIV and 
CIII] broad lines. In addition, they identified the H+K break 
at ~4100A and Mgb absorption (5381 A). This spectrum was 
not of sufficient quality to identify any further features, but 
was su ccessful in i t s prim ary aim: to confirm the system as 
a lens. iFoltz et al.l (|l992l ) observed the galaxy to determine 
a velocity dispersion for the bulge. They identified Mgb and 
Fel absorption features, but also only presented a spectrum 



over a small wavelength range (5100-5550A). iLewis et al.l 
(1998) concentrated their observations around the quasar 
emission lines to investigate variations between the four im- 
ages, thereby inferring microlensing effects (line and contin- 
uum regions are expected to behave differently under mi- 
crolensing conditions because the flux emanates from re- 
gions of different source size). The system 2237+0305 has 
also been the subject of several microlensing monitoring ex- 
periments to study t h e size of the quasar emitting regions 



perimcnts to study tnc size 01 tne quasar emitting regions 
(Woz niak et al. 2000: ; lUdalski et al.1 Booll lEigenbrod et all 



2008a). Rauch et al.l lj2002l ) identified absorption lines in the 



spectrum from intervening (between lens and source) objects 
using the HIRES spectrograph on Keck. These spectra dis- 
play high-resolution absorption line profiles of the interven- 
ing systems and as such do not concentrate on the spectrum 
of t he galaxy as a whole- 
Ivan de Ven et al.l (|2008l ) have recently obtained CMOS 
IFU data of the Call near-infrared triplet from the central 
region of 2237+0305 to determine the contribution of dark 
matter in the central regions of the galaxy. They demon- 
strated that within the central projected 4", dark matter 
can contribute no more than 20% of dynamical mass, us- 
ing single stellar population models to constrain the stellar 
mass contribution, and both lensing and inferred kinemat- 
ics to constra in the total dynamical m ass. Our results agree 
with those of Ivan de Ven et al.l l|2008l ) in the region of over- 
lap. However, our methodology and aims are significantly 
different, as our goal is to construct a simply parametrized 
model to describe the data over a tenfold radial extent. 

In a previous analysis of the system, ITrott fc Websterl 

(2002) combined lensing and photometric data to constrain 
the inner mass distribution of the galaxy. In addition, neu- 
tral hydrogen rotation points at large radius were used to 
study the influence of the outer dark matter halo. They 
found that additional kinematic information was required 
in the image region to break the remaining model degenera- 
cies. Specifically, the inner logarithmic slope o f the dark mat- 
ter halo, a contentio us issue (see for example lHavashi et al.l 
120041 ; Ide Blokl 12005), could not be constrained without ad- 
ditional information about the inner regions of the galaxy. 

In contrast, this paper presents mass-to- light ratios for 
the bulge and disc components that do not require stellar 
population-synthesis models to constrain the stellar M/L 
ratios, and are therefore independent measurements. The 
only assumption is that the mass-to-light ratio is constant 
across the component. This also allows us to quantitatively 
assess the contribution and density profile of the inner dark- 
matter halo in this high-mass spiral galaxy, complement- 
ing work do ne in dwarf and low surface brightness spi- 
ral galaxies (Ide Blok et all l200ll; iBorriello fc Saluccil l200ll ; 



Ide Blokll2005l; IZackrisson et al.ll20o3 r Furthermore, whereas 
in Trott fc Webster ( 20021 ) the halo was constrained to have 



a softened isothermal shape due to the absence of kinematic 
data, in this paper we have allowed for a more flexible model. 

The paper is organized as follows. Section [2] describes 
spectroscopic observations of 2237+0305 with ESI on Keck, 
presents the resulting optical spectrum and uses absorption 
features to derive both an optical rotation curve and line- 
of-sight velocity dispersion profile along perpendicular axes. 
Section introduces the modelling: Section fc.ll describes the 
mass density models used for each component; a two-integral 
axisymmetric model of the galaxy that is used to generate 
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Table 1. Observing 



the observed kinematics is introduced in Section 13.21 then 
Section ^. 3l describes the algorithm used to find the solution. 
Sections [4] and [5] present the results of the modelling and 
some discussion, and Section[6]concludes. We assume f2 m = 
0.3, tt A = 0.7 and H = 70 kms^Mpc^ 1 throughout this 
paper. 



2 SPECTROSCOPY AND STELLAR 
KINEMATICS 

In this Section, the Echelle Spectrograph and Imager (ESI) 
on the Keck-II telescope is used to measure a high spectral 
and spatial resolution spectrum of the galaxy 2237+0305 
and the lensed quasar images. The spectra are used to ex- 
tract a rotation curve and velocity dispersion profile for 
the galaxy, based upon measurements of the Mgb absorp- 
tion line. The results are consistent with measurements 
l|van de Ven et ai1l2008l ) in the region of overlap but extend 
them considerably in extent. 



2.1 Observations 

The lens galaxy 2237+0305 was observed as part 
of the Structure and Dynamics Survey (LSD, e.g. 
iKoopmans fc TrevJl20o3 ). a project that aims to obtain stel- 
lar kinematic information in lens galaxies to improve mod- 
elling and break degeneracies in their mass models. The data 
were taken on 2001 July 21 and 23 with the Keck II telescope 
on Mauna Kea, Hawaii. Conditions were clear and the seeing 
was 0'.'6. The Echelle Spectrograph and Imager (ESI) instru- 
ment was used in echelle mode, providing a constant pixel 
size of 11.9 kms _1 pixel _1 across all orders, and a spatial 
resolution along the slit of approximately 0.154 arcsec/pixel 
(depending on the order). The galaxy was observed with a 
1.25" slit width giving a kinematic resolution of Av=36.1 
kms -1 . 

The target was observed along two perpendicular slit di- 
rections, sampling close to the lens galaxy major and minor 
axes and trying to avoid the locations of the lensed quasar 
images. Table [T] summarises the key observational data for 
the galaxy. Note that the major axis is at a position angle 
of 77° (measured as usual East from North). Although the 
observations were taken at 12° from the major and minor 
axes to minimise quasar contamination, for simplicity the 
two slit angles will be referred to as the 'major' and 'minor' 
axes herein. The misalignments are taken into account in 
the modelling. Two of the spectra were observed with the 
galaxy offset 5" from the slit centre in order to sample the 
outer regions on one side. 

In addition to the galactic exposures, template stellar 



spectra are required to measure the rotation and velocity 
dispersion of the galactic spectral features. The bulge of the 
galaxy is primarily observed in these observations, and so 
template spectra are needed of old red stars expected in such 
a n env i ronme nt. The I-band images of 2237+0305 presented 
in I Yeei (|l988h show the central 12" to be dominated by old 
and r ed stars with (#-r)=0.58 (V - R=0.9. 1 Windhorst et ail 
Il99lf) consistent with Giant K0III-K5III stars. A set of ap- 
propriate stellar templates, described bv lSand et all i|2004l ). 
were observed during the same run and are used in this 
analysis. 



2.2 Optical spectrum and lens galaxy kinematics 

Data reduction was performed using the EASI2D soft- 
ware package, de veloped by David J. Sand and T. Treu 



l|Sand et al.l 200J). Details of the data reduction are given 
in Appendix |A1 

Figure Q] displays the first eight of the ten spectral or- 
ders (the remaining two orders do not show any relevant 
features of the lens galaxy) obtained by averaging the cen- 
tral four arcseconds in order to increase the signal-to-noise 
ratio. Four arcseconds correspond approximately to the re- 
gion over which S/N > 5 per unbinned spatial pixel. A 
relative flux calibration was performed. Prominent spectral 
features are marked on the spectra. The line-o f-sight ab- 
sorption systems reported in iRauch et al.l l|2002h , between 
the lens galaxy and background quasar, are also confirmed 
in these spectra (z=0.5656: Mgl 4466A, Mgll 4377A, Fell 
4070A; z=0.827: Mgll 5108A, Fell 4867A). There are no 
obvious galactic emission lines at the redshift of the lens, 
as expected from an old stellar population. The Ha absorp- 
tion line (682lA observed) is partly obscured by the B-band 
atmospheric water vapour absorption line in lower resolu- 
tion spectra, however our higher resolution observations can 
separate the two lines and Ha emission is not found to be 
present. The broad emission lines of the lensed quasar im- 
ages, CIV (4175A), CIII] (5145A) and Mgll (7541 A) are 
prominent, wi th CIV and Mgll s uffering from sky absorp- 
tion (see also iHuchra et al.lfl985T ). 

The spectra were used to measure the rotation curve 
and line-of-sight velocity dispersion profile for the major and 
minor axes as described in Appendix [X] Figure [2] displays 
the best fits and their uncertainties. The major-axis velocity 
dispersion profile is relatively flat, immediatel y suggesting 
an u nderlying isothermal density profile (e.g. iBertin et al.l 
l2002h . 



3 THE GALAXY MODEL 



In lTrott fc Websterl (2002), it was argued that kinematic in- 
formation could break the remaining degeneracies (see dis- 
cussions above) between models. No kinematic information 
for the galaxy centre was available at that time, however, 
and degeneracies remained in the modelling. With new high 
quality kinematic information in hand, we can now further 
constrain the mass distribution. 
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Figure 1. Sky-subtracted spectra of 2237+0305 taken at a position angle 12° clockwise from that of the major axis. The spectra have 
been flux calibrated up to a constant factor based on ESI standard response curve, measured in July 2001. Two orders are shown in each 
plot (odd and even orders in red and blue respectively). Relevant features are labelled for the lensed source (s) and for the lens galaxy 
(1). A set of bad columns in order 2 and the atmospheric A and B band in orders 7 and 8 are indicated by shaded regions. 



3.1 Mass models 

Given the larger number of constraints with the addition 
of kinematic information, one can accom modate more flexi- 
ble m odels for the mass components than lTrott fc Websteii 
(2002). There is no reason to change the modelling of the disc 
(exponential surface density) and bar (Ferrers ellipse — al- 
though the non-axisymmetric component of the bar will not 
be considered for the inner kinematic model), since these are 
based on visible light profiles. For the kinematic analysis, the 
exponential dis c is deprojected to an exp onential term for 
the z-direction (|Binnev fc Tremainell 19871 ). The distribution 
of dark matter remains unkno wn, and a more general p rofile 
can be fitted than was used in lTrott fc Websterl (|2002h . 



3.1.1 Stellar component 

For tractability, we assume a model for the bulge component 
that is simple and retains the key features that are required 
both for the lensing and kine matic constrain ts. In previous 
work, and following that of ISchmidtl (| 19961) . the bulge is 
modelled by a flattened de Vaucouleurs surface brig htness 
profi le. The latter is one of a class of Sersic profiles (|Sersid 
ll96S() : 

J( J R) = / exphfe m ( J R/r i) ) 1/ " 1 ], (1) 

rt is the characteristic scale length and b m is a constant. 
A de Vaucouleurs profile corresponds to m = 4, with an 
appropriately scaled rt,. 
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Figure 2. Velocity dispersion profile and rotation curve for the 
major (red squares) and minor (blue triangles) axes. The larger 
error bars far from the galactic centre reflect the lower signal-to- 
noise in these regions. 



iLima Neto et alj (|l999h and lLokas et al.1 (120051 ) present 
de-projections of the circularly symmetric Sersic profile, 
which allow one to use kinematic expressions for the ro- 
tation speed and line-of-sight velocity dispersion. We gen- 
eralise these for a non-spherical density profile, because a 
flattened bulge component is important for the lensing. The 
deprojected elliptical Sersic profile for general m is given by, 



P(0 = Po(ip XP [-(i) 1/m 
T(2m) 

P0 - E °2ar[(3-pM' 
p = 1.0-0.6097/m + 0.05463/m 2 
a = r b /(b m ) m 



(2) 

(3) 

(4) 
(5) 



where a is a generalised scale length, Eo is the central sur- 
face mass density (Eo = ftf,E cr for connection with the 
lensing, and E cr is the critical surface mass density), and 
£ 2 — r 2 + z 2 b 2 /a 2 is the elliptical radius, where a and b 
ar e the semi-ma j or and -minor axes. We use the expressions 
of iNoordermeerl l|2008l ) to describe the rotation curve of an 
elliptical Sersic profile. 

The scale length and ellipticity of the bulge were 
measured from a two-dimensional fit of the HST NIC- 
MOS image (Proposal ID: 7495, PI: Falco). After sim- 
ulatin g the NICMOS p oint spread function with TINY- 
TIM i|Hasan fc Burrows! 1 19951 ). a Sersic plus point source 
(quasars) model was fitted with the galfit software 
(|Peng et al.ll2002T l. yielding a Sersic index of m=4.0±0.1, 
a scale length of rf,=4.07±0.10" and a projected axis ra- 
tio of (a/fc)i). P r =0-64±0.01. We use an intrinsic axis ratio of 
(a/6){,,i„=0.6 (|van de Ven et alj|2008h . 



3.1.2 Dark matter component 

The dark matter halo is modelled with a spherical gen- 
eralised cusp, with variable inner and outer density slope 
jZhadll996l : iMufioz et ai]|200ll ; lK~eetonll200ll ). This profile 
reduces to a variety of useful models, such as the NFW, 
Hernquist and isothermal, and has the form, 



p(r) = 



(6) 



(r/r fe )T[l + (r/r fe ) 2 ](»-T)/2 ; 

where 7 and n are the inner and outer logarithmic slopes, 
respectively, and th is the characteristic scale length. For a 
spherical model, the lensing convergence, k,, is given by, 

'n-1 I s 



k(x) = K h B I - 



§ ) (1+^) (1 - 



n)/2 



2 Fx 



(7) 



2 '2'2'l+a; 2 _ 

where Kh = p s rh/E cr , x = r/rh, B{. . .) is t he beta functio n 
and 2-Fi(. • •) is a hypergeometric function (|Bucholj|l969l ). 
In our model, the outer logarithmic slope, n, will not be 
well constrained as most of the constraints lie within 5" of 
the galaxy centre. It is set to a value of 3, consistent with 
N-body simulations l|Navarro et al.lll996l : lMoore et al.lll999l : 
iNavarro et"alll2004h . 



3.2 Kinematic modelling 

To model the lensing and kinematic data-sets of 2237+0305, 
we assume that the central regions of the galaxy can be well- 
approximated by an axisymmetric two-integral distribution. 

We solve the two-integral Jean's equations from which 
we obtain the velocity dispersion profile and mean stream- 
ing motion of the galaxy. W e follo w the prescription of 
Ivan der Marel fc van Dokkuml (|2007T ). details of which are 
described in Appendix[Bj and solve the two-integral axisym- 
metric Jean's equations on a grid. We model the bulge+bar 
as one single axisymmetric component, ignoring the second- 
order non-axisymmetric bar component. 

Jean's equations provide the first two dynamical mo- 
ments, 



p(x, y, w)v%!(x, y, w)dw, 



E(x,y) 

where w denotes the line-of-sight direction and 



= COS ( 



(8) 



(9) 



> sin 2 i 4- v% 



are the first and second moments, i is the inclination angle, 
and R and cj> are the radial and azimuthal directions. The 
rotation component, V, and velocity dispersion, a, are then 
derived as, 



V 

2 



< v w >, 



V 2 



2 

" °~RMS 

<vl> 



(10) 



(11) 



The effects of seeing, pixel size and slit width are in- 
corpora ted into the modellin g using the Monte-Carlo tech- 
nique of Ivan der Marelll| 19941) . The luminosity- weighted mo- 
ments are convolved with the point spread functions that 
characterise the slit width, pixel size and seeing. This is ac- 
complished by randomly sampling the region of the galaxy 
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Component 


Parameter 


Value 


Definition 


Bulge 


K b 


free 


convergence 




n> 


4.07" 


scale length 






64 




Rotation 


k 


free 


rotational support 


Disc 


h '<j 


free 


convergence 




I'd 


11.3" 


scale length 




a/b d 


0.5 


axis ratio 


DM Halo 


Kh 


free 


convergence 






free 


scale length 




7 


free 


inner slope 




n 


3.0 


outer slope 


Bar 




free 


convergence 


Source 


Px 


free 


x source position 




Py 


free 


y source position 



Component 



Parameter 



Value 



Table 2. Parameters to be fitted in the combined lensing and 
kinematic model of 2237+0305. The ellipticities of the bulge and 
disc are fixed at their photom etric valu e , and the bar and disc 
scale lengths are modelled as in ISchmidd jl996h with only a vari- 
able M/L. There are 9 parameters in total. 



in question and averaging luminosity-weighted samples: a 
point is chosen randomly (uniform sampling) from within 
the slit and pixel, and then a normal deviate (with FWHM 
equal to the seeing) is added. The kinematics at this point 
are computed and added (with a weight proportional to the 
galaxy luminosity at that point) to previous samples. This 
method produces good results without having to perform ex- 
pensive numerical convolutions. We sample 3000 points for 
each datum to arrive at our observed moments. By comput- 
ing kinematics with a larger number of samples and com- 
paring output, we find that 3000 samples are sufficient to 
produce accurate results (error <10% of measurement er- 
ror). 



3.3 Mass- model optimization 

The algorithm for consistently solving for the kinematics 
and lensing constraints involves a simplex minimisation of 
the x 2 °f the mass model. Specifically, we: (1) construct a 
galaxy mass model as described in Section T3. II with parame- 
ters shown in Table [2] (2) calculate the lens image positions 
based on this model; (3) calculate the observed rotation and 
velocity dispersion profile based on a summation of the disc, 
bulge and dark matter halo components; (4) calculate the 
expected rotation at the location of the outer HI measure- 
ments from all mass components; (5) compute a \ 2 f° r the 
model as a summation of the adequacy of fits for the lens 
positions, inner kinematics and outer kinematics; (6) itera- 
tively change the parameters until an acceptable minimum 
is obtained. 

The parameters in the model are given in Table [21 
including the unknown source position. W e use t he HST- 
derived image positio ns of I Crane et all (|l99ll ). as m 
iTrott fc Websterl (|2002l 1. 

The constraint^] include the four image positions, the 



1 Flux ratios were not used to constrain the models due to the 
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6.6±O.3T 
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=4.07" 
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Disc 


M/L B 

Td 


1.2±O.3T 
0.014±0.004 
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DM Halo 


Kh 

rh 
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Table 3. Parameter values for the best-fitting solution with re- 
duced x 2 = 2.2. Here, Pae and Pan refer to the source position 
relative to the galactic centre. Uncertainties are also shown cor- 
responding to an increase of one in the x 2 , while marginalising 
over the remaining parameters. 



two HI rotation points and the kinematic information, with 
points in the minor axis rotation curve removed due to con- 
tamination, giving 54 constraints, and 45 degrees of freedom. 

The x 2 statistic quantifying the goodness-of-fit for a 
model is defined as, 
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where the summations are over lens image positions, HI ro- 
tation points, Keck major and minor axis rotational values 
and Keck major and minor axis velocity dispersion values, 
respectively. The reduced x 2 i given the large number of de- 
grees of freedom, is defined as x 2 /45. A reduced x 2 of unity is 
a statistically acceptable fit. Image position errors are 8mas, 
consistent with measurement errors. 



4 RESULTS 

Table [3] lists the parameters for the best- fitting model and 
the 1 a uncertainties for each parameter, defined as the 
parameter value corresponding to an increase in the x 2 of 
one while marginalising over the other parameters. This is 



added complexities of modelling microlensing, substructure and 
the effects of differential dust extinction, however we do calculate 
the resultant ratios and compare them with mid-IR measurements 
as a model sanity check. 
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Figure 3. (Upper) Rotation curve (edge-on) for the solution with 
the best \ 2 . Light-shaded regions indicate radii where kinematic 
information is available. The dark-shaded region in the galactic 
centre indicates where the lensing information is useful. The two 
HI points are shown with their la measurement uncertainties. 
The contribution of the luminous component alone (disc+bulgc) 
is also displayed. (Lower) Rotation curve for a solution with min- 
imum dark matter halo. This solution is allowed within the 3ct er- 
ror bars and demonstrates the disc-halo degeneracy in the model. 



achieved by sampling select points in parameter space, and 
determining where the \ 2 increased by one. Specifically, to 
determine the uncertainty on parameter, alpha, for exam- 
ple, we performed a full minimization for several values of 
alpha (i.e., set the value of alpha and allow all other param- 
eters to vary to find the minimum) and find the value for 
alpha where the \ 2 value increased by one from the global 
minimum. 

Figure |3] (upper) displays the circular velocity curve for 
the best model. Also displayed are regions of the galaxy 
where kinematic or lensing information adds constraints. As 
expected, the model with the best-fitting solution is bulge- 
dominated in the central regions and has a stellar disc that, 
combined with the bulge, supports the majority of the ro- 
tation tor~ 20", beyond which the dark matter halo dom- 
inates. The inner dark matter halo slope is consistent both 
with an NFW profile and a central core at 3cr. The bulge 
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Figure 4. Degeneracy between the disc mass and halo virial 
mass, with contours at 1, 2 and 3<r above the minimum (mini- 
mum x 2 +2-3, 6.2, 9.2). 



and disc mass-to-light ratios are consistent with previous 
measurements. Table [3] also displays the image magnifica- 
tions and total magnification for the model. The flux ratios 
are co nsistent with the mid-IR measurements of lAgol et al.l 
(2000), providing a good sanity check for the model. 

Figure[3] (lower) displays the rotation curve of a solution 
that is at the very edge of what is allowed by the uncertain- 
ties (3cr), having a minimum dark matter halo and a more 
prominent disc. The relatively large errors on the HI data 
points, the freedom of the dark matter model and the lack 
of kinematic information between 5 and 30 arcseconds, com- 
bine to allow at 1% level a low-mass halo model. This demon- 
strates that the familiar disc-halo degeneracy, where the disc 
convergence and shear can account for the halo contribution, 
is almost completely broken. Figure [3] demonstrates the de- 
gree of degeneracy between the disc mass and halo virial 
mass. As expected, these parameters are correlated, but the 
system does not allow a large range for either. In a system 
with a larger disc-halo degeneracy, the range of disc and 
halo kappa values would be increased. In 2237+0305, the 
bulge contributes the majority of the lensing convergence 
and shear, with the disc and bar both also contributing con- 
vergence and shear, and the halo contributing convergence 
alone (spherical model). Although the bulge cannot change 
to account for the disc shear, there is freedom in the bar 
model to account for this. As such, 2237+0305 is a good 
system to break the disc-halo degeneracy. Even tighter con- 
straints can be obtained in the future if kinematic informa- 
tion is available at a larger range of radii (notice that the 
shape of the total rotation curve is different between the 
two figures). Future studies, including a prior on the stel- 
lar mass-to-light ratio of the bulge and improved kinematic 
data, should be able to further limit the remaining degener- 
acy between the halo and disc mass. 

We emphasize, however, that there are correlations be- 
tween other modelled parameters. The bulge, being the dom- 
inant component for the lensing, is well-constrained and dis- 
plays only mild correlations with both the disc and halo 
kappa (results not shown), supporting our assertion that 
the disc-halo degeneracy is the most prominent. In terms of 
halo parameters, the halo convergence (kh) and inner slope 
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Figure 5. (Left) Major (upper) and minor (lower) axis data (pink squares) and fits (blue diamonds) for the rotation curve (inclined 
and calculated over an aperture) for the best-fitting solution. The measurement uncertainties on the data are also shown as error bars. 
(Right) Major (upper) and minor (lower) axis model and data line-of-sight velocity dispersion profiles for the best-fitting solution. 



(7) display an anti-correlation, as do and scale length 
(th)- Although these correlations do exist, the small error 
bars on the bulge convergence tightly constrain the overall 
model. We note that there is freedom in the halo model, but 
this is to be expected given that there are no independent 
constraints on its shape. This freedom in the halo model is 
reflected in the relatively large error bars quoted on its pa- 
rameters. By marginalizing over remaining parameters, the 
error bars include the effects of these correlations between 
parameters. 

The model fits to the kinematics are also of interest. 
Figure [5] shows the fits to the inner regions. The major axis 
kinematics are well recovered, within the observational er- 
rors, but the minor axis dispersion and rotation appear to 
be systematically low. It is possible that the omission of the 
bar in the central kinematics has contributed to this discrep- 
ancy. 

Figure [6] displays the critical lines and caustic for the 
best-fitting model. The swallow-tail caustic feature and con- 
sequent non-elliptical critical line are due to the highly ellip- 
tical bar component. Note that the source does not lie close 
to the swallow-tail caustic and thus the details of the bar 
component do not strongly affect the image positions and 
magnifications. 

Overall, the best-fitting model reproduces the data well 
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Figure 6. Critical lines and caustic for the best-fitting model. 
Also shown are the image positions (asterisks) and source po- 
sition (plus). Note that both radial and tangential caustics are 
displayed. 

considering the coarseness of the data and models. The 
kinematics are reproduced well with a kinematic best \ 2 
of 76 for 44 kinematic data points (with 7 parameters to 
be fitted, giving 37 degrees of freedom). The lensing re- 
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Component 


Mass enclosed 




xlO 8 M 


Halo 


8.8±2.6 


Bulge 


142.5±7.0 


Disc 


4.1±1.0 


Bar 


11.6±1.9 



Table 4. Mass enclosed within the Einstein radius (0.88") for 
each of the four mass components. 

duced x 2 is somewhat higher than unity, as expected be- 
cause in general simple smooth models are not expected to 
reproduce the image configuration with astrometric preci- 
sion. However, this level of precision is more than sufficient 
to determine ro bustly the total enclosed mass within th e 
Einstein radius (|Kochaneklll99ll ; iTreu fc Koopmansll2004h . 
the physical quantity that we are mostly interested in for 
our model. The mass within the Einstein radius (0.88") is 
(1.67=b0 .07) x 10 10 h7oMo, consistent with previous measure- 
ment s (fWambsgan ss fc PaczvnskH 1 19941 ; Ivan de Ven et all 
2008). The components of this enclosed mass are described 
in Table [4] The components of the \ 2 are as follows: the 
lens image positions contribute 21.0% (representing 14.8% 
of the constraints), the HI points contribute 3.0% (3.7%), 
the inner rotation contributes 20.6% (42.6%), and the inner 
dispersion contributes 55.4% (38.9%). 

4.1 Disc maximality and global parameters 

With the best fitting model in hand we can now try to ad- 
dress the question of whether the stellar disc of the lens 
galaxy 2237+0305 is maximal or not. The relevant quan- 
tity for comparison with other observations and theory is 
V2.2/Vtot, i.e. the ratio between the disc rotation veloc- 
ity and the total rotation velocity measured at 2.2 scale 
lengths, close to the peak of the disc rotation curve. Based 
on the best fitting model, and taking into account the er- 
rors on the best fitting parameters, we find that V2.2/Vtot 
= 45±11%, consistent with a sub-maximal disc. The best 
fitting model also allows us to derive the global proper- 
ties of the galaxy. The virial mass of the dark matter 
halo (calculated in a radius with an average density of 
200/Ocrit) is (2.0±0.6) x 1O 12 M0 (concentration parameter, 
c = r 2 oo/r h = 13.3), the stellar mass is (7.4±0.3) x 10 10 M Q 
for the bulge and (7.6±1.9) x 10 10 M Q for the disc. As ex- 
pected for this particular lensing configuration, the bulge 
component is very well constrained, but interesting results 
can be obtained on the disc as well. The stellar to total mass 
ratio is 7.0 ± 2.3 % in good agreement with typical numbers 
for galaxies of comparable mass. 



5 DISCUSSION AND COMPARISON WITH 
PREVIOUS WORK 

The constraint s used to model the galaxy in 
iTrott fc Websterl (|2002l ) were image positions and the 
two neutral hydrogen rotation points at large radius (and 
well away from the image positions). Given that no kine- 
matic information was available in the central regions, the 
inferred velocity curve due to the four mass components was 



assumed to be circular. This means that an acceptable fit 
to the data was attainable with the four smoothly modelled 
mass components because the constraints effectively lie at 
only two radii — 0.9" (images) and ~30" (HI data). There 
were considerable dege neracies in the model at other radii 
(e.g. see Figure 2 from lTrott fc Webster|[2004l ). 

In the present analysis, the addition of kinematic data 
provides additional constraints out to a radius of ~5" (~4 
kpc) . As displayed in Figure [3] constraining information is 
available in the central 5" and at two outer points in the 
galaxy, leaving the models to fit the remainder of radii ac- 
cording to their input density profiles. The bulge is well- 
constrained because it contributes considerably to the mass 
within the inner 5". The mass-to-light ra tios of this study 
are m ore tightly constrained than those of lTrott fc Websterl 
(2002) — their la uncertainties were calculated while hold- 
ing the values of other parameters at their minimum, 
thereby underestimating the true uncertainties. The bulge 
mass-to-light ratio is consistent with previous measure- 
ments from early-type galaxies M/L_B=7.9±2 .3Tp),fc = o.7o. 
ilGerhard et all l200ll; ITreu fc Koopmanj 12004m . Similarly, 
I Jiang fc Kochanekl l|2007h report M/Lb = 7.2±0.5T Q , h=0 .70 
from an analysis of 22 lens galaxies. As far as the disc mass- 
to-light ratio is concerned, our best fitting value is also con- 
siste nt with typical measur ements reported in the literature 
(e.g. ICourteau fc Rixlll999l . and references therein). 

The dark matter halo, however, contributes mostly in 
the outer regions of the galaxy where few constraints are 
available. In addition, unlike the bulge and disc, the ab- 
sence of a light profile to guide the density profile of the 
dark matter yields additional freedom in the model. This 
allows a more flexible model, but also increases degeneracy. 
It is thus remarkable to see how well the dark matter halo 
is constrained by the available data. However, more detailed 
information about the structure of the dark matter halo 
would require information at radii greater than 6", where 
its influence increases. Alternatively, to gain more informa- 
tion about the halo from these data, one could place a prior 
on the bulge and disc mass-to-light ratios based on studies 
of other galaxies (or on stellar synthesis models) and deter- 
mine whether such constraints narrowed the allowed value 
of the inner slope. This is left for future work. 

The best-fitting solution described in lTrott fc Websterl 

(2002) had a more massive disc and less massive halo com- 
pared with the best-fitting solution described here. The 
bulge mass was consistent between the two studies, within 
quoted errors. The errors quoted in lTrott fc Websterl (2002) 
were calculated assuming no parameter covariances, and 
therefore underestimate the true errors as calculated by 
marginalising over the remaining parameters (as performed 
in this work). Therefore it is difficult to determine whether 
the halo and disc masses (and therefore the disc mass-to- 
light ratio) are consistent between the two studies. Note, 
however, that the dark matter halo mass profile used in this 
study was different to that in ITrott fc Websterl (|2002l ). al- 
lowing freedom in the inner slope. 

The fit to t he image positions in t his study is not as 
good as the fit in lTrott fc Websterl (|2002l ). This is due to the 
increased complexity of the model in this work, and the need 
fo r the fit to balance the im age positions and the kinematics. 
In lTrott fc Websterl (|2002l ). the image positions formed the 
most precise portion of the data (the HI points, being far 
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from the galactic centre and with relatively large error bars 
are easier to fit), and therefore were able to be well- fitted. 
In this work, the first and second order kinematics formed a 
relatively larger proportion of the \ 2 calculation, and there- 
fore were fitted at the expense of the image positions. 

In future work the modelling could be improved to in- 
clude any effects of external shear (G2237+0305 has other 
less massive systems identified along the line-of-sight) and 
to include the kinematic effects of the bar more consistently. 
The former could improve the lensing component of the \ 2 
without affecting the dynamical fit. In addition, the dark 
matter halo could be modelled to be non-spherical with flat- 
tening aligned with that of the disc (as suggested by studies 
of other lensing systems) . 

Despite the comparatively good model presented in this 
present work, it is important to keep in mind the following 
caveats. The lens galaxy 2237+0305 has additional structure 
that is not accounted for in the four component model, e.g. 
prominent spiral arms that are attached to the bar compo- 
nent. In addition, the bar is not treated in a self-consistent 
manner. These may not be important for the lensing but 
still produce anisotropy in the kinematics. The bar is dif- 
ficult to model kinematically and an attempt has not been 
made in this work. However, in the central few arcseconds 
its presence is small and we opted to neglect it for the inner 
kinematics, including its effect only in the lensing and for 
the overall rotation curve. 

The image magnifications and total magnification are 
dependent on the balance between the steepness and the 
anisotropy of the overall mass m odel. As demonstrated by 
IWambseanss fc Paczvnskl J 19941 ). 2237+0305 allows a range 
of models with degeneracy between the density slope and el- 
lipticity. Inclusion of kinematics in our model has reduced 
that degeneracy. The magnificati on ratios are con sistent 
with the mid-IR measurements of lAgol et all (|200Qt) . They 
argued that mid-IR measurements are the least likely to 
be contaminated by ste llar microlensing, electron scattering 
and e xtinction (see also lChiba et al . 2005; lAgol fc Kochanekl 

l20oeh . 

In terms of global properties, the central galaxy is bulge- 
dominated, and modelling the mass distribution of the bulge 
accurately is most important for an accurate overall model. 
Beyond r ~20", the halo appears to continue to increase in 
rotational support in the outer regions, and to absorb the 
rotation beyond the optical radius. 

Our results can be compared with those of 
Ivan de Ven et all (|2008f ) . Their measured central ve- 
locity dispersion (166 ± 2 kms -1 ) is consistent with that 
presented in this work (172 ± 9 kms~ ) within the errors. 
Although they had the advantage of two-dimensional data 
and the additional information that accompanies that, 
our data are based on stellar kinematics from the Mgb-Fe 
region, observed in a cleaner region of the sky than the 
Call triplet targeted by their observations. This allowed 
us to reach further out in projected distance from the 
centre, albeit along pr eferred directions. The results of 
Ivan de Ven et ail (|2008h are complementary to ours: they 
fit a fourier-based lens model, fitting the image positions 
and radio flux ratios, and combine the information with 
the implied mass distribution from the kinematics and 
light distribution, to obtain an estimate of the luminous 
mass-to-light ratio. A comparison of this value with that 



obtained using a single stellar population model yields a 
maximum dark matter contribution within the Einstein 
radius of 20%, which is higher than our best fit value of 
7% and our upper limit of 15% (3<r). In addition, their 
comparison of the lens-based mass model and the observed 
luminosity distribution indicates that mass follows light in 
the inner regions, and any halo contribution is relatively 
constant. Our analysis focuses on determining a simple yet 
comprehensive description of the galaxy over a broader 
range of scales. Although the results for the dark matter 
halo can be improved by more extended and precise kine- 
matic data, we achieved our aim of determining the bulge 
and disc mass-to-light ratios with good precision, without 
the need to use stellar population synthesis models. Future 
studies with the addition of stellar population synthesis 
models will be very valuable. 

The good agre e ment between our results and those 
of Ivan de Ven et al.l l|200Sl ) - based on largely independent 
datasets (except for the strong lensing) and methodologies 
- is encouraging: the combination of lensing and dynamics 
is a powerful tool to study spiral galaxies, not only elliptical 
galaxies. As larger samples o f spiral lens galax ies are be- 
ing discovered by the SLACS (|Bolton et al.1 12008'1 and other 
surveys - including many where the Einstein radius will be 
equal or larger than the scale length of the disc - the com- 
bination of lensing, kinematics, and stellar population syn- 
thesis models should enable substantial progress in our un- 
derstanding of the structure of disc galaxies, and therefore 
their formation and evolution. 



6 SUMMARY 

Data from the ESI echelle instrument on the Keck II tele- 
scope were used to measure the major and minor axis kine- 
matics of the spiral lens galaxy 2237+0305. The rotation 
and velocity dispersion profile in the central few arcseconds 
of 2237+0305 are combined with the lensed-image positions 
and the rotation curve beyond the optical disc (HI rota- 
tion points) to construct a structural and kinematic model 
of 2237+0305. Four components are used to constrain the 
mass distribution — a dark matter halo with variable inner 
logarithmic slope and scale length; a de Vaucouleurs bulge, 
modelled structurally with a flattening consistent with the 
light profile; an exponential disc, inclined at 60°; a Ferrers 
bar, at a position angle 39° from the disc major axis. The 
inner galaxy kinematics were modelled with an axisymmet- 
ric two-integral model composed of the bulge, disc, and dark 
matter halo. 

The best-fitting solution shows that the lens galaxy is 
bulge-dominated in the inner regions, with a massive stellar 
disc and dark matter halo that provide rotational support 
beyond the bulge. The bulge and disc B-band mass-to-light 
ratios (M/L 6 =6.6±0.3T©, M/L d =1.2±O.3T ) are consis- 
tent with typical values found for local galaxies. The bulge 
mass is very well constrained by the combination of lensing 
and dynamics. The disc mass is less well constrained be- 
cause of residual degeneracies with the halo mass, although 
the best solutions indicate a sub-maximal disc, contributing 
45+11 % of the support at 2.2 scale lengths. The inner log- 
arithmic slope of the dark matter halo is consistent with a 
negative inner logarithmic slope of 7 = 1, with a preferred 
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value of 7=0.9±0.3 (note that this may still be dependent on 
our use of an axisymmetric two- integral kinematic model) . 

This work illustrates the potential of lensing and dy- 
namics t o investigate the interna l structure of spiral galaxies 
(see also Ivan de Ven et al.1120081 ). The system Q2237+0305 
was the natural starting point for this kind of investigation 
due to the unusually low redshift of the lens galaxy. It is 
encouraging that such a complex galaxy, including a bar in 
the central regions, can be described with a relatively simple 
model described here over a wide range of scales. A partial 
limitation of this target, however, is that the multiple images 
are located at a particularly small physical scale - signifi- 
cantly smaller than the effective radius of the bulge and the 
exponential scale length of the disc - due once again to the 
unusually low redshift of the lens galaxy. This implies that 
the lensing geometry has relatively little leverage on the dark 
matter halo, compared to typical lenses at higher redshifts 
where the Einstein radius can be up to a few half-light radii, 
where the dark matter halo starts to be a substantia l frac- 
tion of the total mass (e.g. iTreu fc Koopmansl lix)! ) . Cur- 
rent and future samples of disc lens galaxies are likely to be 
found at higher redshift, thus providing an opportunity to 
break the disc-halo degeneracy even further, as well as to 
investigate evolutionary trends. 
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APPENDIX A: DATA REDUCTION 

The ESI observations require some special preparation and 
calibration before they can be combined and used for kine- 
matic analysis. The calibration process undertaken includes 
the following steps: bias subtraction, flat-fielding, rectifica- 
tion, cosmic-ray rejection and sky subtraction. These steps 
were performed by the package EASI2D, which is d eveloped 
by David J. Sand and T. Treu (|Sand et al.1 120041 ) for easy 
extraction of echelle orders. 

Given the lack of emission lines in the spectra of the 
old stellar population, strong absorption features need to 
be used in the kinematic analysis. The canonical region for 
stellar absorption kinematics in the optical - the Mgb-Fe 
complex around 5200 A rest frame shown in Figure IA1I - 
was used as our primary kinematic dataset. As an additional 



test, we also fitted the NaD doublet at 5892A, which is gen- 
erally consi dered less reliable b ecause of possible interstellar 
absorption l|Sparks et al.ll 19971 ) . 

In practice, the Gauss-Hermite Pixel Fit- 
ting Software developed by R. P. van der Marel 
( |http: / /www-int.stsci.edu/^marel /software/pixfi t.html[ | is 
employed to determine the kinematic properties along both 
the major and minor axes. The template and galaxy spectra 
are prepared to have the same resolution, number of pixels 
and wavelength range. They are the n compared us i ng th e 
iterative fitting method described in Ivan der Marell l|l994h . 
The variable parameters used to find the best-fitting (low- 
est x 2 ) solution are the order of the polynomial fit to the 
continuum and the spectral type used as the template. The 
rotation curves were found to be quite robust to parameter 
variations, while as expected the velocity dispersion profile 
requires more care in choosing the parameters. The region 
shown in Figure IA1I was found to be a good compromise 
between using a large enough region and avoiding as much 
as possible the broad AGN features. A large number of 
tests were run to explore spectral ranges, masking areas 
and polynomial order for continuum fitting. Our error bars 
include all these sources of errors which dominate over the 
random errors. 

The radius range is also extended beyond the region 
where the galaxy signal is obvious in order to include as 
much information as can be obtained from the data. The 
data are binned radially to 0.6" bins, in order for the points 
to be independent given the seeing. The outer points do 
not produce good fits — beyond r~4", the galaxy signal 
is lost in the noise. The best-fitting stellar templates are 
from Giant stars with typ e G9II I - K5III, consistent with 
the imaging results of I Yed (|l988l ) and produce statistically 
indistinguishable fits. Figure [AT] displays an example fit. The 
stability and accuracy of the kinematic results are very good 
— all templates produce mutually consistent fits and the 
residuals appear consistent with noise. 

The NaD region gave consistent rotation curves to the 
Mgb-Fe region, but systematically higher velocity disper- 
sion profiles by 25%. Possible explanations of this difference 
- in addition to the effects of interstellar absorption - are 
contamination by quasar light and the diversity of possibly 
kinematically distinct stellar populations in the targeted re- 
gions. Although we feel that our choice of the Mgb-Fe region 
is justified, we repeated our analysis using the NaD-based 
stellar velocity dispersion profile and found our results to 
change only marginally. None of the conclusions of this pa- 
per is affected significantly by this choice. 

The centre of the galaxy is difficult to locate from the 
data due to the significant contamination from the lensed 
quasars and the difficulties associated with image deconvo- 
lution. We therefore symmetrise the data to match the shape 
of the curve on both sides to find the kinematic centre. Once 
the centre is found, the rotation curves are symmetric within 
errors. The kinematic profiles obtained with the four best 
fitting templates have been averaged to produce the final 
major axis results. The analysis of the minor axis data is 
similar to that performed for the major axis. 

The central line-of-sight velocity dispersion is consistent 
within la uncertainties between the major and minor axis: 
Oc ^ 172±9 kms" 1 , whic h compares well with the result 
of Ivan de Ven et all (|2008l ) of a c ~ 166±2 kms" 1 . Our un- 
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Figure Al. Best kinematic fit in the spectral region of the Mg 
and Fe absorption lines (rest wavelengths are indicated on the top 
axis of the plot). The data are shown as a black histogram, while 
the best fitting template is shown as a red curve. Spectral regions 
affected by intervening absorption line systems are masked out 
during the fit and shown as shaded areas. 



certainties are larger, due to our inclusion of dominant sys- 
tematic errors such as template and continuum mismatch. 
Our results are m arginally lower than the value inferred by 
iFoltz et all (|l992h (215 ± 30 kms" 1 ), possibly for the same 
reasons that cause the mismatch between the Mgb-Fe and 
NaD results. 



APPENDIX B: TWO-INTEGRAL 
AXISYMMETRIC GALAXY MODELS 



Here 



briefly 



summarise 



the 

20071 ) to 



method of 
produce ax- 
See 



Ivan der Marel fe van Dokkuml 

isymmetric dynamical models from photometric data 
the original paper for further details. 

We assume the galaxy can be fitted with an axisym- 
metric model, and that there are two integrals of motion. In 
this regime, the two coupled Jean's equations are: 



dpv 



9$ „ 

oz oz 



dpv 2 



9$ 



OR 



where v 2 R 



0. 



(Bl) 
(B2) 



The gradients of the gravitational potential, V3>, are 
calculated from the currently used mass density profiles 
of th e mass components, according to (|Binnev fe Tremaind 
1 19871 ): 



N 

F = -V$ = -vrG^yi 

=1 



Pi(m )Vm 



2/2 R 

m /a,i — 



r + af 



+ 



+ bV 



(B4) 



and a,i,bi are the semi-major and -minor axes of the mass 
component, and the summation is over all distinct mass 
components. 

The mean strea ming motion is not constrained by the 
Jean's equations, and Ivan der Marel fe van Dokkuml (|2007l ) 
introduce a constrained scaling parameter, k, to control the 
amount of rotation, 



(B5) 



The value of k is constrained by the necessity of positive 
dispersion everywhere, 



k < minK/(u?-u|)] 



1/2 



(R,z) 



(B6) 



The mass components and velocity moments then need 
to be inclined to the viewing angle, and the projected line- 
of-sight moments determined according to, 



{v$){x,y) 



1 



p(x, y, w)v$(x, y, w)dw, 



E(s,y) 

where w denotes the line-of-sight direction, and, 



(B7) 



(B8) 



v w = Vtf) cos <p sin z, 
Vw = — v %) cos 2 <f> sin 2 i + v'ft 

are the first and second moments, and i is the inclination 
angle. The rotation component, V, and velocity dispersion, 
a. are then derived as, 



V 

2 



= < V w >, 
If 2 2 

— y —crms 



= V 2 



(B9) 



(B10) 



dr, (B3) 



where, 



The prescription to determine the first and second ve- 
locity moments using this method is as follows: 

• For a given mass model, calculate the gravitational po- 
tential derivatives and density profiles on a (R, z) grid; 

• Calculate v 2 = on the same grid according to IB2I 

• Calculate v 2 ^ on the same grid according to IB2I 

• Calculate the mean streaming motion on the same grid, 
according to lB51 and ensure it meets the positivity require- 
ment; 

• Incline the density distributions of the mass compo- 
nents to the line-of-sight; 

• Calculate the line-of-sight first and second velocity mo- 
ments; 

• Convolve the moments with seeing, slit width and pixel 
size. 
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